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Abstract. The observed power spectrum of high-degree solar p-modes [t > 200) shows discrepancies with the 
power spectrum predicted by the stochastic excitement and damping theory. In an attempt to explain these 
discrepancies, the present paper is concerned with the influence of the observed subsurface flows on the trapped 
acoustic modes (p-modes) . The effect of these inhomogeneous background flows is investigated by means of a non- 
modal analysis and a multi-layer model. It is shown that the rotational and meridional components of the velocity 
field change the wavelengths of the oscillation modes which, in turn, results in modifications of the corresponding 
modal frequencies. The magnitudes of the frequency residuals depend on the spatial scales of the modes and on 
the gradients of the different components of the flow velocity. Together with other mechanisms (e.g. the scattering 
of modes by the large scale convection (Goldreich & Murray 19943), the non-modal effect of the variation of the 
frequencies in time may contribute: 1) to the observed widening of the corresponding peaks in the observed power 
spectrum with increasing angular degree; 2) to the partial dissipation of spectral power, and, as a result, 3) to 
the discrepancies between the predicted and the observed power spectrum of solar p-modes. 
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1. Introduction 

One of the main goals of helioseismology is to understand 
the excitation and damping mechanisms that yield the ob- 
served power spectrum of the solar oscillation modes. The 
theory of the random generation of modes by convective 
turbulence has been developed by Goldreich et al. I| 1994(1 

n 



and Goldreich & Kumar (199U|) in an attempt to quan- 
titatively understand the structure of mode sources and 
the generation rate of the oscillation modes. Observations 
show (see e.g. Libbrecht ({T^55j) : Woodard et al. (|2001[l ) 
that the energy of the p-modes with frequencies below the 
photospheric cut-off grows with the angular degree £ for 
intermediate values of £, i.e. 100 < I < 200. This obser- 
vational fact is in good agreement with the predictions 
from the theory of turbulent excitement and damping of 
p-modes mentioned above. However, for £ > 200 there 
are substantial discrepancies between the theory and the 
observations. As a matter of fact, recent observations by 
Woodard et al. ( 2001) show that for I > 200 the wave en- 
ergy distribution function decreases rapidly with increas- 
ing angular degree I. The latter authors remark that the 
discovered decline in mode energy at high i values may 
indicate the existence of some 'unmodeled mechanism of 
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damping'. In the present paper, the viability of such a 
mechanism is studied, viz. mode damping due to a sheared 
background flow. As a matter of fact, there exist numer- 
ous helioseismological measurements of solar subsurface 
flows (see Gonzalez Hernandez et al. (|1999fl . Gonzalez 
Hernandez & Patron J. l^OUOfl h These observations show 
that the velocity field of these flows has both rotational 
and meridional components and that the gradients of these 
flow velocity components depend on the spatial coordi- 
nates. We will present a multi-layer model for these flows 
that allows to quantify their effect on the solar p-modes. 

Several other candidate mechanisms for non-turbulent 
wave damping have been suggested, e.g. the nonlinear cou- 
pling between trapped and propagating modes (Kumar & 
Goldreich (1989)) or mode scattering by convective mo- 
tions (Goldreich & Murrav ri994f) . Moreover, the direct in- 
fluence of the inhomogeneous temperature, the magnetic 
field, and/or the velocity fields on the solar p-modes and /- 
modes have been studied intensively (Murawski & Roberts 
l)1998al[T998b|l : Murawski & Goossens ifTMH)) : Vanlommel 
& Cadez (fTMSf . Vanlommel & Goossens JEHU), while 
'resonant absorption' (or 'continuum damping') has been 
introduced as a candidate damping mechanism for p- 
modes by Tirry et al. (1998). Clearly, the physical char- 
acteristics of the solar medium affect and systematically 
distort the main frequencies and the travel times of the 
modes. This insight has led to the development of the 
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time-distance technique in helioseismology (Duval et al. 
(1993)) which became a powerful tool for the observation 
of 3D structures with inhomogeneous magnetic and tem- 
perature profiles. 

As pointed out by Goldreich et al. I|1994[) . when the 
widths of the frequency peaks in the power spectrum are 
greater than the inverse of the observation time, these 
peaks contain information about the damping rates of the 
modes. These authors also indicated that the excess width 
of the peaks (with respect to the inverse observation time) 
can result from a variation of either the frequency or the 
amplitude of the modes and that it is reasonable to as- 
sume that both these types of velocity modulation yield 
comparable contributions in the line width. Observations 
show that the line widths of the oscillation modes increase 
with increasing frequency and angular degree. Goldreich 
& Murray 1)1994[) suggested the scattering of modes by 
turbulent velocity fluctuations as a possible mechanism 
supporting this observational fact. 

It is well known that the main properties of the so- 
lar oscillation modes are well described by the standard 
normal mode analysis. For a static equilibrium and for 
appropriate boundary conditions, the linear MHD (mag- 
netohydrodynamics) operator is self-adjoint. However, in 
neutral fluids and in plasmas with nonuniform background 
flows there are often discrepancies between the analyt- 
ical results of the standard modal analysis and exper- 
imental data. These discrepancies are usually ascribed 
to the mathematical incompleteness of the normal mode 
spectrum due to the non-self-adjoint nature of the gov- 
erning equations (Trefethen et al. I|1993|) . Criminale & 
Drazin (2000), (1990)). One of the alternative and com- 
plementary methods to solve the initial value problem 
in sheared plasma flows is the non-modal analysis. In 
this formalism a transformation of the variables is per- 
formed and, as a result, the spatially inhomogeneous 
terms in the governing equations are replaced by time- 
dependent terms. Therefore, in general, this method en- 
ables to study the non-exponential temporal evolution 
of linear perturbations of sheared flows. The non-modal 
analysis was introduced originally by Lord Kelvin ( 1887) 
and the theoretical basis for studying the wave dynam- 
ics in hydrodynamical and plasma flows has been de- 
veloped by Goldreich & Linden-bell (|1965(l . Craik & 
Criminale (1986); Chagelishvili et at. I|1997|l : Rogava et 
al. (2000] I2001[) . Many authors showed that the shear 
in the flow velocity field can have substantial effects 
on the wave dynamics, e.g. causing time-dependent fre- 
quencies and wave numbers (Butler & Farrell (Q992); 
Chagelishvili et al. H1994fl ) and coupling of different wave 
modes (Chagelishvili et al. H199hJ) : Rogava & Mahajan 
(1997)). The non-modal approach has also been applied 
to the study of several aspects of the excitement of gravi- 
tational and acoustic waves on the Sun (see Chagelishvili 
et al. I|2000|) : Pataraya & Pataraya l|1997[> . 

The significant discrepancies between the observed and 
predicted power spectra of trapped modes with high an- 
gular degree I > 200 may be due to sheared flow effects. 



As a matter of fact, these modes mostly propagate in the 
thin cavity below the photosphere (r > 0.95i?Q) in which 
plasma flows with strongly pronounced velocity gradients 
have been observed. Hence, the aim of our study was to in- 
vestigate the possible role of the non-exponential temporal 
evolution of the modes (due to the nonuniform flows) in 
the formation of the observed power spectra. We applied 
a non-modal analysis to study the effect of the observed 
sub-photospheric shear flows on the frequencies and the 
propagation characteristics of the solar p-modes. In order 
to distinguish the effect of velocity shear from all other 
effects, we adopt a very simple model and consider the 
modes locally, in plane slab geometry. Our model allows 
to separate the effects related to the spatial inhomogene- 
ity of the background flow from the inhomogeneity coming 
from the gravitational stratification and the temperature 
gradient and enables us to focuss on the evolution of the 
wave characteristics associated with the flow inhomogene- 
ity itself. 

In the next section, we present the model we ap- 
plied to investigate the interaction between the inhomoge- 
neous flows beneath the photosphere and trapped acous- 
tic modes. In the third section, we describe and analyze 
our theoretical results. In the fourth section, we discuss 
the obtained numerical results and we present and discuss 
some phenomenological points of view in the framework 
of the obtained results. Then, we discuss the possible ob- 
servational consequences of the shear flow effects on the 
p-modes. In the final section, we formulate our conclu- 
sions. 

2. Physical model and equations 

The aim is to study the temporal evolution of small per- 
turbations in a stratified medium with a non-uniform two- 
dimensional 'background' flow. The problem is considered 
in plane-parallel geometry and is described in Cartesian 
coordinates. The ir-axis denotes the 'toroidal' direction 
and its positive direction is chosen as the direction in 
which the Sun rotates. The y-axis coincides with the 
'meridional' direction and is directed towards the north 
pole. The z— axis covers the 'radial' direction and is di- 
rected outwardly, to the solar atmosphere. We study the 
properties of modes with small amplitudes which are con- 
sidered as linear perturbations of the equilibrium state. 

2.1. The equilibrium model 

In this subsection, the main parameters characterizing the 
equilibrium are introduced, in particular the parameters 
characterizing the equilibrium flow velocity. These param- 
eters are determined in such a way that they yield a good 
approximation of the above-mentioned measurements of 
the solar subsurface flow profiles. In general, the subsur- 
face flows can be very complicated. In the non-modal anal- 
ysis, the nine components describing the shear rates of the 
background flow velocity Vo components in the three spa- 
tial directions arc collected in the so-called 'Shear Matrix ' 
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Fig. 1. Schematic view of the modelled subsurface spher- 
ical shell (in plain slab geometry). The dashed horizontal 
lines schematically represent the boundaries of the thin 
layers considered in our multi-layer model. Within each 
of these thin layers the sound speed is assumed to be ap- 
proximately constant and they correspond to the piecewise 
constant temperature profiles shown in Figure or [3 



In order to consider the propagation of small perturba- 
tions in the gravitationally stratified medium we adopt a 
local analysis, similar to the formalism used by Whittaker 
(1963) and Ulrich Ijl97()[> . It should be noted, however, 
that we here focuss on the effects of the inhomogeneous 
background flow and thus we excluded the effects related 
to the magnetic field, radiative flux and convective mo- 
tions. The latter issue has already been addressed by sev- 
eral authors (e.g. see Swisdak and Zweibel H1999|) V In ad- 
dition, the influence of the random velocity fields on /- 
modes has been studied by Murawski & Roberts (|l993a 
1993b|) and by Murawski & Goosscns (1993) also includ- 
ing the effect of the chromospheric magnetic field. Within 
each of the thin shells in our model (schematically shown 
by the horizontal dashed lines in Figure 1 the tempera- 
ture (i.e. sound speed) is assumed to be constant. Hence, 
the temperature profile of the standard solar model (see 
e.g. lStandard Solar Model 455|) is approximated by a step- 
function. In other words, locally, in each thin layer we 
have an equilibrium density and pressure profile of the 
form pq(z),pq{z) oc exp (— z/if), where H is the local scale 
height, so that in each layer the temperature is constant. 



S (for details e.g. see Mahajan & Rogava 
most general form this matrix reads: 

V()x : x V()x,y VQx,z 

S = | V0y iX Voy^y Vfty <Z 

Vqz,x V$z,y Vo z . z 



). In its 



(1) 



where the index after the comma denotes the derivative 
with respect to the corresponding spatial coordinate. We 
want to study the joint influence of the rotational and 
the meridional subsurface flows on the solar p-modes. 
Hence, the equilibrium velocity field includes both rota- 
tional (x— ) and meridional (y— ) components. In our gen- 
eral model we consider the x-component of the equilib- 
rium flow velocity as a function of both the y- and z- 
coordinates, Vq x — Vo x (y, z). Further, we also take into 
consideration the meridional component of the flow veloc- 
ity and its dependence on the radial coordinate, i.e. Voj, = 
Vo y (z). Consequently, in the equilibrium state we have an 
inhomogeneous velocity field Vq — (V~o x (y, z), Vq v (z), 0). 

In section 4 we will model the subsurface flows with a 
linear approximation of the velocity profiles measured by 
Gonzalez Hernandez & Patron H2()()()J) . Hence, in the local 
frame of reference the components of the flow velocity read 
as: 

V 0x = Vbi + ay + bz, (2) 



and 



V Qy = V 2 + cz, 



(3) 



where, Vbi and V02 represent the constant parts of the 
rotational and meridional components of the flow velocity 
and the parameters a = Vo Xiy , b — Vq x%z and c = Vo VyZ are 
the local spatial derivatives of the corresponding velocity 
component. Note, that the considered background flow is 
incompressible V • Vq = 0. 



2.2. The temporal evolution of the p-modes 

Let us now study the propagation of small perturbations 
in the considered plasma shell and focuss on the effects 
related to the inhomogeneity of the background flow(s). 
We consider the waves as linear perturbations of the sta- 
tionary equilibrium described in the previous section. The 
following set of linearized equations then governs the dy- 
namics of the small disturbances in a gravitationally strat- 
ified isentropic gas: 



dpi 
dt 



■ (Vo-V)pi + (v r V) Po (z) + p (z)(V-i>i) = 0, (4a) 



dvi 



Po(z)— + p (z) (Vq ■ V) V! + p (z) ( Vl ■ V) V 

ot 



-Vpi + pig, 



(4b) 



dpi 
dt 



(V -V) Pl + (vv 1 V)po(z) = 



7Po ( dp_ 
dt 



Po 



- + (Vo-V)pi + {v v V)p {z^j , 



(4c) 

where g denotes the gravitational acceleration and 7 is 
the ratio of specific heats, symbols with index denote 
the equilibrium quantities and those with index 1 the per- 
turbations. 

Clearly, in the case of a static medium, i.e. without 
the background flows, we get a simplified set of equations 
without the inhomogeneous terms coming from the con- 
vective derivatives. In this case, the normal-mode solution 
of this set of equations takes the form: 



if> ~ exp 



i I Lot + k x x + k y y + k 



-) 

2H J 



(5) 
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where, ip represents the quantities pi/po, Pi/pa and the 
components of the velocity perturbation. However, we 
want to know what happens to these solutions when the 
non- uniform equilibrium flows are taken into account. 
Notice that Eqs. (0} contain convective derivatives, which 
also take into account the spatial inhomogeneity of the ini- 
tial (equilibrium) velocity field. These spatially inhomo- 
geneous terms complicate the analysis considerably. For 
small magnitudes of the background flow and of its gradi- 
ents the problem has been studied in the standard frame- 
work of the normal modal analysis and the influence of 
the flow on the eigenfrequencies has been treated as a 
small perturbation of the basic frequencies. This pertur- 
bation technique is applied in Backus and Gilbert (1961) 
and Ulrich et al. I|1979|l . The methods of inversion of the 
solar internal rotation rate, from the frequency splitting 
observations of the solar non-radial oscillations, are based 
on a similar theoretical formalism. These methods allow 
to determine the internal rotational velocity profile of the 
Sun with a high accuracy. 

However, due to the inhomogeneous background flow 
the governing equations have a non-self-adjoint character. 
Hence, we apply the non-modal technique, which enables 
us to find an additional class of 'non-normal' solutions of 
the above set of equations 10} . These solutions describe the 
evolution of the perturbations in time which, in general, 
can be non-exponential. 

Following the non-modal technique we represent the 
perturbation quantities pi/po, Pi/po, and Vi in the fol- 
lowing form: 



ip{x, y, z; t) = if>(k x , k y , k z ; t) exp [i \ip l - ip 2 - — 



2H 



(6) 



where, 



and 



ipi(x,y,z;t) = k x (t)x + k y (t)y + k z (t)z, (7) 



Mkx(t),ky(t);t) = V i J k x (t)dt + V Q2 J k v (t)dt. (8) 

In these expressions k x (t) , k y (t) and k z (t) are the (in gen- 
eral) time dependent components of the wavevector satis- 
fying the following set of ODEs: 



dk 
~dl 



S T -k = Q, 



(9) 



where S T is the transposed shear matrix. In this general 
formalism, the substitution of expressions ©-© in the 
governing equations results in a set of ordinary differen- 
tial equations (ODEs) in time in which the spatial inho- 
mogeneity is replaced by a temporal one (for details see 
e.g. Rogava et al. ( 2001)). These ODEs govern the tempo- 
ral evolution of the spatial Fourier harmonics with time- 
dependent components of the wave vector. In the case of 
the two dimensional non-uniform flow considered here, one 
gets: 

f=0, (10a) 



dky 
~~dT 



ak x = 0, 



dk z 
dt 



bk x + ck v = 0. 



(10b) 



(10c) 



Consequently, for the considered flow velocity profiles the 
IE-component of the wave vector remains constant in time, 
kx = k x0 , but 

k y (t) = k yQ ~ ak x t. (11) 

Integrating Eq. i|10c|) . we find for the temporal evolution 
of the radial (z-) component of wavevector: 



k z (t) = k z0 - (bk x0 + ckyo) t + ackx0 t 2 . 



(12) 



It should also be noticed that (Rogava 122221) there 
exists a combination of the solutions of the Eqs. (|10a|l - 
(|10cf) that is preserved in time, viz. 



A = ck 2 . + 2(bk y k x — ak z k x ) = const. 



(13) 



Using the representation @ in Eqs. ||2J) and rewriting the 
equations for the components of the velocity perturba- 
tions, normalized density perturbation D — pi/ po and the 
enthalpy perturbation Q = pi/po) we get the following set 
of ODEs for the perturbed quantities: 



dD 
~~dt 



+ i [k x v x + ky(t)u y + k z (t)v. 



-Jp + av y + bv z = ~ik x Q, 
^ + cv z = -ik y (t)Q, 



dQ = C 2 f dD , ( 7 - l)v z 
dt s \ dt 



0. 



(14a) 

(14b) 
(14c) 
(14d) 

(14e) 



where k z (t) = k z (t) + i/2H. Hence, the presence of shear 
flows causes a transformation of the wave solutions in 
time. The rate of this transformation strongly depends on 
the initial orientation of the wave vector and on the values 
of the shear parameters included in the shear matrix. 

Our first aim is to qualitatively understand the tempo- 
ral behavior of the p-modes driven by the joint effect of the 
inhomogeneous rotation and the meridional flow appear- 
ing in the thin shell just below the photosphere. Therefore, 
in the next section we will introduce some simplifying as- 
sumptions. This will enable us to derive an approximate 
analytical 'dispersion relation' which governs the tempo- 
ral evolution of the 'effective frequencies' of the modes. 
We will also demonstrate the shear flow effect in a simple 
model with a linear velocity profile, viz. for a mode with 
Z = 951 which is confined within a very thin cavity above 
the fractional radius 0.997. In section 4 we then try to 
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quantify this effect in first order with a multi-layer model. 
The measured velocity profiles are approximated linearly 
in each layer and the cumulative effect of the layers is cal- 
culated taking into account the time spent by the modes 
in each layer. 

3. Approximative analysis 

In general, the application of the non-modal approach 
means that, in the governing equations, the spatial in- 
homogeneity (related to the terms containing the nonuni- 
form background flows) is replaced by a temporal one. 
After this transformation of the perturbation variables the 
temporal behavior of the spatial Fourier harmonics (con- 
taining the time dependent wavenumbers) can be studied. 
In most cases, the solution of the resulting ODEs (in time) 
requires the application of numerical techniques. In this 
particular case of the trapped acoustic-gravity modes, the 
initial value problem has a specific feature. On one hand, 
the time dependence of the wavevector results in a vari- 
ation of the phase speed and the direction of the wave 
propagation, in time. On the other hand, the temperature 
gradient refracts the propagating mode. Formally speak- 
ing, the local approach outlined in the previous section 
is the 'non-modal equivalent' of the so called ray approx- 
imation, commonly used for the local description of p- 
modes (Whitaker flMty : Ulrich (JTHZOI))- But, in this case 
we do not perform a Fourier transform with respect to 
time. Instead, one has to study the propagation of the 
modes by solving the governing equations (|14|l in each 
layer separately and then matching the solutions at the 
interfaces between the different layers. However, at this 
stage we aim to qualitatively estimate the characteristic 
changes in the modal properties in order to understand 
the principle of the mechanism of these changes induced 
by the inhomogeneity of the background flows. Therefore, 
we adopt an approximative analysis to study the propaga- 
tion properties of the modes and to construct the shape of 
their ray paths. The part of the ray path confined within a 
given layer, with constant temperature, is locally linear (a 
schematic view of this kind of piecewise straight ray path 
is shown in Figure El and within this layer the equations 
governing the motion of a point of the wave front can be 
written as follows: 



upper t.p. 



upper t.p. 



Sx \ , . kx 



Tt)r vphM w\ 



m(t)k Xi (zj,t) 

\h(z t ,t)\2 ■ 

Lj(t)k Zi (Zj,t) 

\ki(z h t)\ 2 ' 



(15a) 



(15b) 



(15c) 



where, index i denotes a given layer between horizontal 
planes located at zu — Zi and Z2i = z< + (Sz)f, ki, k Xi , 
k Vi and k Zi are the module and components of wavevec- 
tor, respectively. In these equations we introduce the time 
dependent 'effective' frequency w(t). We will turn to the 




lower t.p. 



Fig. 2. Schematic view of the piecewise straight ray path. 



issue of validity of such a description in the following sub- 
section. It should be noticed that, because of the temper- 
ature gradient the local values of the sound speed in each 
neighboring layer are different. Therefore, the components 
of the wavevector should be adjusted accordingly, in or- 
der to keep the 'effective' frequency unchanged through 
the interface between the layers. One can easily see that 
in the case without a shear flow the 'effective' frequencies 
should be replaced by the usual constant modal frequen- 
cies. In this latter case the problem always can be reduced 
to two dimensions (as it is usually done under the normal 
mode formalism). In our case, however, the trajectories 
(ray paths) of the modes are, in fact, 3D curves. From 
this point our problem is reduced to the determination 
of the 'effective' frequency and how it relates to the time 
dependent components of the wavevector. 

3.1. The 'effective' frequency 

Generally speaking, the effect of an inhomogeneous flow on 
waves manifests itself in a temporal variation of the modal 
properties. Using different analytical methods (among 
which the non-modal approach) it has been shown that 
strong gradients in the background flow velocity can cause 
very rapid changes of the modal frequency and the am- 
plitude, even in time intervals as short as an oscillation 
period. As a result, the power spectrum of the waves is 
affected in such a way that it may become impossible to 
distinguish the individual peaks in the spectrum. 

In the Sun, however, such strong gradients of the flow 
velocity are not observed. Therefore, in this first attempt 
to model this 'small' effect qualitatively, we assume that 
the sheared flows below the surface of the Sun just slightly 
transform the oscillation modes and that this transforma- 
tion is represented by a small deviation of the mode prop- 
erties found by the standard modal analysis. As a result 
of the flow inhomogeneity there appear terms of the form 
SijUi (where a summation over repeated indices is meant) 
in the equations (I14bl - (j!4cfl . Here the indices i, j — 1, 2, 3 
represent the x—, y—, and z— direction, respectively, and 
S is the shear matrix. From Table 1 it is clear that the 
gradients of the flow velocity components are two to four 
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orders smaller than the characteristic angular frequencies 
of the oscillation modes. Hence, we expect that the con- 
tribution from the SijUi terms in the solution of the gov- 
erning equations is very small. Therefore, at this stage we 
neglect these terms, at least in equations (|14bl) - (|14cl) . 

The above approximation results in a mathematical 
simplification of the problem as it makes the frequency 
purely real, i.e. oscillatory. In other words, it is equivalent 
to the assumption that the amplitude of the mode varies 
much slower in time than the rapid oscillation itself so 
that the amplitude remains approximately constant. This 
approximation thus results in perturbations of the form 
i\) w %p(to)exp(i(p(t)). The set of Eqs. i(T4"|) then leads the 
dispersion equation for dip(t)/dt, which has a similar form 
as the one given by Whitaker ( 1963) for waves propagating 
in a static medium: 

(^r) 4 - (^r) 2 W M + ® + ^)+^ - °. 

(16) 

where, k\ = k^+ky, cu a = r yg/2C s is the local acoustic cut- 
off frequency, and uji, = 'g(j — l)/-/H is the local Brunt- 
Vaisala frequency. The difference between Eq. 1)16(1 and 
the dispersion equation for 'homogeneous' (in the sense of 
flows) waves is that now the coefficients in the equation are 
time dependent. Thus by using Eq. © we can write the 
approximate 'effective frequency' as the time derivative of 
the time dependent phase: 



Voikx 



V 02 ky 



dk v 

'-at 



dk z d<p(t) 



at 



dt 



(k ■ V ) - 



d<p(t) 
dt 



where, 



d<p(t) 
dt 



i w . 



-^W^-A^kUt), 



(17) 



(18) 



W 



C^(kl(t) + k z (t) 2 )+col (19) 

In Eq. H18|) we chose the '— ' sign before the square root 
without loss of generality. It can easily be shown that the 
first four terms in the RHS of Eq. (|17fl in fact represent the 
expression (fc • Vo). It is well-known that this scalar prod- 
uct of the wavevector and the background flow velocity 
represents the effect of the Doppler shift of the frequencies 
due to the fact that the waves are advected by the flow. 
Hence, this term is analogous to the one given for example 
by Ulrich et al. (1979). However, here the components of 
the wavevector (and thus the characteristic frequency) are 
time dependent. This time dependency represents the de- 
viation from a purely exponential (~ exp(iwi)) temporal 
evolution of the small disturbances. 

3.2. Effect of a simple (linear) velocity profile 

In order to numerically calculate the changes of the 
modal frequencies purely due to the non-modal effects 
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Fig. 5. Results of the calculations for the same mode as 
in Fig. 0] for a smaller value of the shear parameter a = 
—0.1 • 1CP 6 s _1 . The value of the shear parameter b has 
been taken the same as in Fig. 0] (b = 2.5 • 10~ 5 s _1 ). In 
panel A we show the curve of the frequency variation in 
time, for the angle j3 = 45°, within the first three evolution 
hours. The solid line corresponds to c = 2.0 • 10~ 5 s -1 
and the dashed line to c = 3.0 • 10 _5 s _1 . Similar curves 
are plotted in panel B for the angle (3 — 30°, where the 
dashed-dotted line corresponds to the smaller value of the 
parameter c and the dotted line to the larger one. Panel C 
shows the resulting modal frequency evolution over a time 
interval of 10 hours. The different line styles correspond 
to the line styles in the previous panels. 



we developed a numerical code making use of a spatial 
discretisation on the basis of the Standard solar model 
(www . ap . stmarys . ca/ ~guenther/solar/ssm455 . sink) 
and the approximate 'dispersion relation' (Eq. l|17fh. 
In our calculations we used also the observed 
p-mode frequency data (published online by 
Libbrecht et al. and Bachmann ct al., see the 
web page http://www.gong.noao.edu/teams/data/ 
jwl_freqs.html). As an example we examined the mode 
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Fig. 3. Piecewise constant radial profile of the temperature (panel A) and scale height (panel B) corresponding to the 
case of the simple velocity profile addressed in subsection 3.2. 
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Fig. 4. The variation of the modal frequency in time. The curves correspond to the mode £ = 951, n = 2 and with 
the basic modal frequency vq = 4.5243 mHz. The calculations have been done for a simple, linear profile of the flow 
velocity with the shear parameters equal to a = — 1 • 10~ 6 s — 1 , b — 2.5 • 10~ 5 s^ 1 and c = 2.0 • 10~ 5 s _1 . Panel A: curve 
a corresponds to k y o > (angle (3 — 45°) and curve b corresponds to k y o < (angle (3 = —45°). The points denoted 
by la and 16 correspond to the frequency of the respective mode close to its upper turning points (i.e. the frequencies 
observable on the surface!) while the inner turning points are denoted by 2a and 26. Panel B: curves of the frequency 
variation corresponding to the upper turning points for different directions of the horizontal wavevector k^ (i.e. for 
different values of the angle f3). The curves are denoted respectively as: a x and a 2 - (3 = ±5° (crossed lines), &i and 
6 2 - P = ±30° (dashed lines), C\ and c 2 - fi = ±45° (solid lines), d\ and d 2 - (3 — ±60° (dotted lines), e\ and e 2 - 
(3 = ±85° (dash-dotted lines). 



with I ~ 951. In Figure we show the piecewise constant 
radial profiles of the temperature (panel A) and scale 
height (panel B). The shown range of the fractional 
radius includes the very thin cavity in which the £ = 951 
mode is confined. 

Notice that these illustrative calculations involve a 
sample mode with I — 951, n — 2 and vo — 4.5243 mHz. 
The curves in Fig. 01 correspond to a simple velocity pro- 
file. In particular, it was assumed that the profiles of both 



the rotational and the meridional flows are linear in the 
entire (narrow) cavity in which the mode is confined. As 
one can see from Eq. (11711 the non-modal change in wave- 
length directly manifests itself in a variation of the ef- 
fective frequency in time. However, the problem we are 
considering here has one specific feature: the modes are 
trapped in the acoustic cavity between the upper and in- 
ner turning points. As a result, if the modal frequency 
increases (decreases), when it propagates from the upper 
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turning point to the inner one, then the frequency de- 
creases (increases) when mode propagates in the backward 
direction. Therefore, any observable changes in modal fre- 
quencies (on the surface), can only arise when the net 
influence of the flow, during the propagation of the mode 
in the inward and outward directions, is nonzero. Hence, 
the nonzero residuals in Fig. 0] result from the asymmetry 
of the background flow. The example shown in Fig. ^cor- 
responds to a case with a simple flow velocity profile: for 
this calculation we assumed linear rotational and merid- 
ional flow profiles in the entire cavity. This simple profile is 
appropriate for the modes with very high angular degree. 
The total residual during the observed modal lifetimes is 
thus the sum of the net frequency changes appearing dur- 
ing each passage from the upper turning point to the in- 
ner one and backward as shown in Fig. 0] (Panel A). The 
curve a (solid line) corresponds to the mode with k x > 
and k y o > and the curve b to the k x > and k y o < 
(dashed line). The points la and 2a correspond to the 
moments when the mode a reaches its upper and inner 
turning points, respectively. Similarly we have points lb 
and 2b for the curve b. The Panel B in Fig. ^illustrates the 
changes in modal frequency for different ratios k x /k y (0), 

1. e. for different directions of the horizontal wavevector in 
the initial momentum. The curves on this panel in fact 
are the lines connecting the residuals in frequencies cor- 
responding to the moments when the considered mode is 
close to its upper turning points. We denote the curves 
for modes with k y o > by the characters with the index 1 
and for the modes with k y o < by characters with index 

2. If we represent the k y o components of the wavevector as 
k y o = kh sin(/3), where (3 is the angle between the kh and 
x-axis, then the curves on Panel B correspond to the val- 
ues of (3 in the following order: curves a\ and (crossed 
lines) correspond to {3 = ±5°, &i and &2 (dashed lines) 
correspond to (3 = ±30°, c\ and C2 (solid lines) corre- 
spond to (3 = ±45°, di and (doted lines) correspond to 
(3 — ±60° and e± and ei (dash-dotted lines) correspond to 
(3 — ±85°. Clearly, the residuals in modal frequency due to 
the non-modal effects are greater for modes propagating 
in an oblique direction with respect to the both rotational 
and meridional flows and the effect is maximal for the an- 
gle (3 — ±45° (i.e. k x — k y o) (the angle corresponding for 
maximum rates of the frequency residuals can differ from 
this, say it can be ±30° or ±60°, depending of values of 
shear rates and wavevector components). 

In Fig. El we show results for a rather small value of the 
parameter a (a = —0.1 • 10~ 6 s _1 ) obtained for the same 
mode as in Fig.^J The curves in panel A correspond to the 
angle (3 = 45°. The solid curve corresponds to the shear 
rate c = 2.0 • 10~ 5 s -1 and the dashed line corresponds 
to c — 3.0 • 10~ 5 s -1 . Similarly, we show the results for 
the angle (3 = 30° in panel B, where the dashed-dotted 
curve corresponds to the smaller value of parameter c and 
the dotted line to the larger one. It should be noticed here 
that from both panels A and B one can easily see that the 
amplitude of the frequency variation along the ray path 
increases as c increases. Finally, we combine the results 



for both angle values (3 — 45° (solid and dashed lines) and 
(3 = 30° (dashed dotted and dotted lines) calculated for a 
10 hour period. In this panel the line styles correspond to 
the same values of the shear parameter c as in the previous 
panels. From this panel we conclude that the rate of the 
overall frequency residuals drops down by a factor of 10 
compared with the case shown in Fig. ^ since now the 
value of parameter a is 10 times smaller. This seems to 
indicate that the frequency residual rates depend nearly 
linearly on the value of this parameter. We will come back 
on this point in the next section. 

The aim goal of this subsection was to outline the prin- 
ciple of the process and to get a rough estimate of the 
possible frequency residuals for different parameter values. 
For these purposes we adopted a simplified velocity profile 
with fixed constant values of the shear parameters within 
the entire cavity (above the fractional radius 0.997), in 
which the examined oscillation mode with £ = 951 is con- 
fined. This simple model yields the following main conclu- 
sions: (1) The joint effect of the non-uniform rotational 
and meridional flows can yield finite differences between 
the normal mode frequencies and the 'effective frequen- 
cies'; (2) These 'residuals' can be of the order of a few tens 
of /iHz. It is thus worthwhile to improve the model in an 
attempt to better 'quantify' the effect of the shear flows. 
This is done in the next section, where we consider a more 
advanced model including observed (not simplified) veloc- 
ity profiles, and perform calculations for different sample 
modes including the one considered in the current subsec- 
tion. 

4. Numerical results for observed velocity profiles 

In this Section we consider cases of modes with differ- 
ent angular degrees. For these purposes, in general, more 
complicated, the observed velocity profiles should be taken 
into account in the model. Modes with £ significantly lower 
than that considered in the previous subsection, penetrate 
deeper in the solar interior and correspondingly we expand 
the width of cavity, where we consider the propagation on 
modes as it is shown in Figures|S]and[7| Further we explain 
the model in detail. 

4.1. Modelling the subsurface flows 

As a basis for the present study we have the results of 
helioseismological measurements of the subsurface veloc- 
ity profiles (see e.g. Gonzalez Hernandez et al. ([1999), 
Gonzalez Hernandez & Patron (2000)). These observa- 
tions show that the gradients of the flow velocity com- 
ponents are functions of position. That is the reason why 
we present our results for a number of magnitudes of the 
shear parameters confined within some ranges of values. 
The model 'shear parameters' included in the shear ma- 
trix are estimated as follows. The photospheric latitudinal 
variation of the rotation rate is usually modelled as: 

fi(0)= a k cos k (6)+ So , (20) 

fc=0,2,4 
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0.9981 - 0.9992 
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-80.4253 


0.996 - 0.998 


-14.2857 






0.998 - 0.999 


21.2500 






0.999 - 1.000 


28.5714 



where 6 denotes the colatitude and clq = 452.0 nHz, 
a 2 = -49.0 nHz, a 4 = -84.0 nHz, and s = -31.7 nHz 
(Snodgrass l|1984|l . Gonzalez Hernandez & Patron |2000)). 
As the parameter a ~ Vo x , y included in the shear ma- 
trixes corresponds to the derivative with respect to the 
y— coordinate of the ^-component of the equilibrium flow 
velocity. Thus, this parameter represents the latitudinal 
variation of the rotational velocity in the considered local 
plane-parallel slab. We performed the calculations assum- 
ing that this parameter is of the order of a ~ —0.01 • 10 -6 , 
-0.1 • 10- 6 , -0.3 • 10~ 6 , -0.6 • 10~ 6 , or -1.0 • 10~ 6 s" 1 . 

To estimate the values of the parameter b = Vq x ^ z , 
which locally represents the radial variation of the ro- 
tational component of the flow velocity, and the radial 
shear c of the meridional (y— ) component of the velocity 
field, we used results of helioseismic inversions given in 
(Gonsalcz Hernandez et al. 1^999), Gonzalez Hernandez 
& Patron (2000)). In Fig. we show the rotational profile 
for latitude 15° taken from Gonzalez Hernandez & Patron 
(2000) (dotted line). We have made a linear approxima- 
tion of the observed profile. The approximate profile is 
also plotted in Fig. [fj] (solid line) . This approximate pro- 
file corresponds to specific values of the shear parameter 

b. The results of our estimates, for different fractional dis- 
tances from the solar center, are shown in the Table 1. 
At depths below the thin sub-photospheric layer, i.e. for 
r < 0.95i?Q, we assume the radial gradient of the rota- 
tional velocity to be negligibly small. Here it should be 
noticed that we only consider oscillation modes with an- 
gular degree I > 95, which have their lower turning points 
far above the solar tachocline. These parameters a, b, and 

c, determine the inhomogeneity of the background flow. 
Together with the spatial scales of the waves themselves 
and with the timescales of the interaction between the 
modes and the flows, they determine the strength of the 
influence of the inhomogeneous flow on the waves. 

4.2. Determination of the timescales of influence 

Another important issue we have to address concerns the 
timescales of the interaction between the flow and the 
oscillation modes. One of the basic temporal parame- 
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Fig. 6. Profile of the rotational velocity in the solar sub- 
surface layers. The solid line shows the sample profile for 
latitude angle 15° given by Gonzalez Hernandez & Patron 
(2000). The dotted curve is the observed rotational profile. 
The solid lines show the modelled linear approximations 
of the observed profile for different fractional radii. 

ters characterizing the solar p-modes concerns their life- 
times. Several observational methods have been used by 
numerous authors to estimate the lifetimes Tl of the so- 
lar p-modes, see e.g. the work of of Chou et al. (2001) 
and Chen et al. (1996). There exists a discrepancy be- 
tween the observational measurements of the modal life- 
times. Nevertheless, from the above-mentioned measure- 
ments one can estimate the range of the lifetimes as 2-10 
hours. 

Another characteristic temporal parameter of an oscil- 
lation mode is its travel time r along the ray path, from 
one reflection in the upper turning point to another one 
(see Fig. 1). As it is well known these travel times are ob- 
servable by applying the methods of local helioseismology. 

It is well-known that the radial gradient of the 
temperature and, therefore, of the sound speed inside 
the Sun leads to the refraction of waves. Hence, in 
order to properly calculate the changes in the modal 
frequencies as the modes propagate along their ray 
path, we radially discretize the plain-parallel slab. We 
thus consider a set of very thin shells so that within 
each of these thin plasma layers the local sound speed, 
the cut-off frequency, and the Brunt- Vaisala frequency 
can be assumed to be constant. Therefore, the spatial 
orientation of the wavevector also approximately remains 
constant when the waves propagate within the given 
layer. In our calculations we use the equilibrium pa- 
rameters taken from the standard solar model (see e.g. 
www. ap . stmarys . ca/~guenther/solar/ ssm455 . sink). 
We denote the part of the total ray path confined within 
the given layer by ds — v p hdt where, v p h = to/k is the 
modal phase velocity and dt is the time interval during 
which the wave front stays within the given layer. This 
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Fig. 7. As in Figure |2l corresponding to the range of the 
fractional radius values r = 0.95 — 1.00. The squares in the 
right bottom corners show the areas plotted in respective 
panels of the Figure |3| 

distance relates with the distance between the top and 
bottom boundaries of the layer dz as follows: 

dz = ds cos a, (21) 

where a is the angle between the vector ds and the verti- 
cal (see the schematic view in Fig. QJ. The cosine of this 
angle can be written in terms of the components of the 
wavevector. In particular, cos a = k z /k where k z and k 
are the radial component and length of the wavevector, 
respectively. Finally according to the above expressions 
we can write: 

k 2 

dt = ——dz. (22) 
uik z 

This last expression determines the time interval during 
which the wave undergoes an influence from the inhomoge- 
neous flow within the given thin layer. Hence, in our model 
we consider the influence of the inhomogeneous flows on 



the oscillation modes as the accumulation of the effects in 
the different layers reached by the given mode as it prop- 
agates from its upper to its lower turning point and back. 
Now, as we have parameterized the effective time of influ- 
ence of the nonuniform flows on the p-modes, the question 
that remains is: "how strong is this influence?" . In other 
words, is the total time interval of influence large enough 
for significant changes in the modal properties or not? 
Below, we concentrate on this issue. 

4.3. Modes with angular degree I < 200 

In this subsection, results of calculations are presented 
based on the model discussed in the previous subsections, 
i.e. the more general case with realistic approximations of 
the velocity profiles. The shear parameters used for these 
calculations depend on depth and are given in Tabic 1. 
In Fig. |21 the results of these calculations are shown for 
different values of the shear parameter a. The panels A, B, 
C, D and E correspond to j3 = ±45° and to the values a = 
-1 • 10~ 8 s-1, a = -1 • 10" 7 s-1, a = -3 • 10~ 7 s-1, a = 
—6 ■ 10 -7 s— 1 and a = — 1 • 10~ 6 s— 1, respectively. Again 
the curves labelled by characters with index 1 correspond 
to modes with k y o > 0, while the labels with index 2 
correspond to modes with k y o < 0. The curves a% and 
a-2. (dashed lines) correspond to the mode with angular 
degree £ = 95, n — 4 and a basic modal frequency vq = 
2.3761 mHz. Clearly, the influence of shear flow is more 
pronounced as time proceeds. It turns out to yield changes 
in the frequency of the order of a few /xHz to, at most, 5 
to 20 /xHz after 10 hours (the maximal life time of the 
mode). The values of the shear parameter a used in the 
calculations correspond approximately to the latitudes 1°, 
10°, 25°, 40° and 60°, respectively. The overall frequency 
residual (after 10 hours) for the considered mode is plotted 
in panel F of Fig. |H1 (dashed line) as a function of the 
parameter a. 

4.4. Modes with angular degree I ^ 200 

From the Eqs. (fTTTfl and lfT?|) - (fTr1f it is clear that the rate 
of change of the modal wavelength in time strongly de- 
pends on the modulus (and consequently on the angular 
degree of modes as k\ oc 1(1 + 1)) and the direction of 
the wavevector, i.e. on its horizontal (kh) and radial (k z ) 
components. In particular, the larger the component kh 
of the wavevector is, the more substantial the changes of 
the modal properties are. Therefore, we can expect that 
the deviation of the modal properties (from those of nor- 
mal modes) increase with angular degree. As a matter of 
fact, our calculations show larger residuals from the ba- 
sic modal frequency for the sample modes (see Fig. |HJ: 
I = 263, n = 6, Vq = 4.2336 mHz (curves b\ and &2 (dash- 
dotted lines)), I — 673, n = 2, Vq = 3.8202 mHz (curves c\ 
and C2 (dotted lines)), and the mode considered in section 
3.2 (curves d\ and di (solid lines)). These plots demon- 
strate the fact that the residuals of the modal frequency 



B.M. Shergelashvili and S. Poedts: The effect of subsurface flows on p-modes 



11 



due to the non-modal effects increase with growth of the 
angular degree. The overall residuals after 10 hours ver- 
sus parameter a are also plotted for these sample modes 
in panel F in the respective order of line styles. This latter 
confirms that the frequency residuals depend on the val- 
ues of the parameter a, and this dependence has a nearly 
linear character. 

In addition, it should be noticed that the modes un- 
dergo a more substantial influence of the flow when they 
are closer to their inner turning points. As a result, modes 
with inner turning points within the area with strongly 
pronounced gradients of the flow velocity (i.e. mostly the 
modes with I ^ 200) , change in time more effectively. This 
fact also contributes to the appearance of more significant 
non-modal effects for modes with higher angular degree. 
Comparing the results shown in Fig.[S]for the sample mode 
with I = 951 (solid lines in panels B and E) with those 
obtained with the simplified velocity shear profiles (see 
panel B in Fig. 0] and panel C in Fig. it is easily seen 
that our preliminary estimates given in section 3.2, yield- 
ing that frequency residuals up to values of a few tens of 
/iHz should be expected, were true. For a = — 1 • 10~ 7 s— 1 
the solid line dl in panel B in Fig. |S| should be compared 
with the solid and dashed lines in panel C of Fig. |5j which 
correspond to the angle (3 — 45°. On the other hand, from 
the panel E of Fig. |H1 one can see that (comparing the 
solid curves dl and d2 with the similar curves in panel B 
of Fig. 0J the characteristic frequency residuals with the 
observed velocity profiles are somewhat smaller than those 
obtained with the constant shear rates in subsection 3.2. 
This is due to the complicated radial profile of the ob- 
served rotational and meridional flows leading to alter- 
nating signs of the shear rates along the radial direction. 

5. Discussion 

Formally speaking, non-uniform flows can influence the 
linear wave modes in fluids and plasmas in two different 
ways. On one hand the flow advects the wave front pro- 
ducing a Doppler shift of the modal frequencies (the scalar 
product on the RHS of the Eq. fTTjl l. On the other hand, 
the shear background flow causes (in general) changes in 
the modal wavelength, which in turn results in a tempo- 
ral variation of the characteristic frequencies and/or am- 
plitudes of the linear disturbances (for the case of modes 
considered here see the last term in RHS of Eq. 1|17|) '). 
Under the standard normal mode formalism it is a pri- 
ori assumed that the perturbations evolve as purely ex- 
ponential (~ exp(iwi)) functions of time. As it is well 
known, the Doppler shift of the modal frequencies can be 
described under the normal mode analysis very well. On 
the other hand, as we have shown the temporal changes in 
wave characteristics related to the time dependence of the 
wavenumbers can be studied as well. The normal mode de- 
composition of the perturbations neglects this latter char- 
acteristic of the temporal evolution of the wave modes. 
For the particular case of p-modes trapped in the solar in- 
terior this assumption is true in most cases. Indeed, there 



are no sharp flow velocity gradients observed in the solar 
interior. Even in narrow shells below the convection zone 
and just beneath the photosphere, where significant gradi- 
ents of the flow velocity are observed, the shear rates still 
are by several orders of magnitude lower than the char- 
acteristic frequencies of the oscillation modes. Therefore, 
we expected that the inhomogeneous flows in the subsur- 
face shell, which is under consideration here, could only 
slightly deviate (non-modally) the temporal behavior of 
the modes from that described under the standard normal 
mode approach. The approximative analysis performed in 
this work leads to the conclusion that these expectations 
were true. Nevertheless, the considered non-modal effects 
still turn out to be 'strong' enough to support different ob- 
servational evidences. In particular, from the global point 
of view the non-modal effects can cause the 'blurring ' of 
the global p-mode power spectra and the mixing of the 
modal power with the noise. On the other hand, these 
effects contribute to a slight 'deformation' of the modal 
ray paths. We address these observational aspects of the 
non-modal effects in details in following subsections. 

5.1. The observational aspects of the non-modal 
effects - global view 

Let us now address the important issue about how the ef- 
fects related to the flow velocity gradients manifest them- 
selves in the observed power spectrum of the solar p- 
modes. It is believed that the solar oscillations are ex- 
cited by the stochastic motions in the upper part of the 
convection zone. From the global point of view, the so- 
lar p-modes are observed as global oscillation modes of 
the Sun by measuring the characteristic frequencies of 
the intensity or the velocity field oscillations. The obser- 
vational data are then concentrated in the power spec- 
trum which enables the identification of different oscilla- 
tion modes by fitting the peaks arising in the spectrum 
with the peaks in the analytically obtained Fourier power 
spectrum. The formation of each peak involves the con- 
tribution of many modes with similar properties. The os- 
cillation modes are characterized by their frequency, their 
radial order n, their angular degree £, and their azimuthal 
order m. As it is well known from the theory of stellar 
non-radial oscillations (Unno et al. I1989|l . in the case of 
a spherically symmetric star the eigenfrequencies are de- 
generate with respect to the azimuthal order. This means 
that it is not possible to distinguish (in the observed power 
spectrum) modes with the same radial order and angular 
degree but propagating in different horizontal directions. 
In Fig. El panel (a) schematically shows one peak corre- 
sponding to a given radial order and angular degree and 
containing the power of all modes with azimuthal order 
m= -£,-£+ 1,...,0,...,£-1,£. 

Any phenomenon that spoils the spherical symmetry 
(for example the axi-symmetric rotation of the star, mag- 
netic field etc.) results in lifting the degeneracy with re- 
spect to m. In the non-spherically symmetric case, instead 
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of one peak their appear several peaks in the power spec- 
trum corresponding to the different orientations of the 
horizontal wave vector. On panel (b) of Fig. lathis is illus- 
trated schematically by two peaks in the power spectrum 
corresponding to two modes with the same radial order 
and angular degree, but with horizontal wave vectors di- 
rected to opposite directions. This effect is the well-known 
phenomenon of 'frequency splitting'. Clearly, frequency 
splitting observations allow to perform the inversion of 
the observational data into the rotational profile inside 
the Sun. It should be noticed, however, that the eigen- 
frequency splitting can also be caused by any other factor 
breaking the spherical symmetry such as, for example, the 
meridional flows. 

The effect of inhomogeneous flows on the power spec- 
trum can be complicated. On one hand, the flow velocity 
gradients contribute to the frequency splitting just like 
the uniform flow would do. On the other hand, we have 
shown above that the influence of the flow inhomogeneity 
can be noticeably important, able to serve as a mecha- 
nism causing changes in the modal properties which can 
contribute to the observed deficit of modal power for very 
high degree p-modes. In the non-modal approach we have 
described these changes by a time dependent wavevector 
and a time dependent modal frequency. Hence, to describe 
the effect of inhomogeneous flows on the power spectrum 
we need to use additional characteristic temporal parame- 
ters: viz. the widths of the peaks in the power spectra Av 
and the observation time, which relates to the resolution 
of the observations. If the changes in the modal frequen- 
cies 5v occurring due to the non-modal effects during the 
lifetime of the mode are much smaller than the half-width 
of the corresponding peak in the power spectrum, then one 
can assume that the changes in frequency are negligible. 
In this case the background flow only contributes to the 
observed frequency splitting. This is the case considered 
under the formalism of the standard modal analysis by 
Ulrich et al. (197HJ). The schematic view of this situation 
is given on panel (c) in Figure 

When 5v ~ Av the changes in frequency contribute 
to the line width (along with different damping mecha- 
nisms causing a variation of the mode amplitudes). When 
Sv > Av the frequencies of the individual modes are 'drift- 
ing' along the frequency axis. In this case a part of their 
power is mixed with the acoustic noise (leading to a de- 
crease of the signal to noise ratio) and a part of their 
power contributes to the neighboring peaks. Under these 
circumstances, even in the conservative case where we as- 
sume that the average power exchanged by two neigh- 
boring peaks are approximately equal to each other (in 
the real situation this can not be true, see Goldreich et 
al. Hiyy^ V the non-modal effects could be the cause of 
the partial dissipation of the power because of the mix- 
ing with noise. The last two situations are schematically 
displayed on panel (d) in Fig. [5] Here we should make 
one important remark. The modes that contribute to a 
given peak in the power spectrum are randomly excited 
by different convective sources at different time intervals. 



But, as the flow velocity gradients are small on the Sun 
the changes in frequencies of the particular mode by the 
non-modal effects need a significant time to occur. That 
is why the peaks in the power spectrum remain to be con- 
centrated around the basic central frequency and why the 
frequency splitting is still observable. 

5.2. The observational aspects of the non-modal 
effects - local view 

Now we turn to another issue, viz., how the non-modal 
changes in modal properties affect the observational data 
obtained by using the well-known methods of local he- 
lioseismology (such as time-distance techniques, ring dia- 
grams, etc.). In fact, the expressions (|15fl determine the 
shape of the trajectory of a point with a constant phase, 
within the framework of our model with a piecewise con- 
stant temperature. One can easily see from these ex- 
pressions that in the shearless (zero shear matrix) limit 
the horizontal components (k x and k y ) remain constant. 
Therefore, in this case the modes propagate in the hori- 
zontal direction along straight lines (x and y coordinates 
of the point are linear functions of time) and the ray path 
is confined in one plane (the problem is two dimensional) . 
The variation of the radial component of the wavenumber 
k z arises from the refraction of waves by the tempera- 
ture gradient along the radius. The schematic view of a 
given modal ray path between two neighboring reflections 
at the upper turning points is shown in Figure 1101 (panel 
A). The effects which can systematically distort the ray 
paths (when we exclude non-modal changes) relate only 
to the inhomogeneous temperature and/or magnetic field 
profiles and Doppler shift of frequencies due to the advec- 
tion of modes by the flows (this effect is represented by 
the scalar product on the RHS of Eq. (|T7|0 . These system- 
atic distortions are immediately detectable by using local 
methods of helioseismology. On the one hand, observa- 
tions detect changes in modal travel times and distances. 
On the other hand, changes of modal frequencies result in 
a transformation of the so-called acoustic ring shape from 
a ring to an elliptical one. 

If we turn now to the case of the nonzero shear rates, 
it can be seen from Eq. (|15fl . that the components of the 
wavevector are time dependent, in general. The charac- 
ter of this dependence is determined by the profiles of 
background non-uniform velocity field. In the particular 
case, which we consider in this paper (rotational veloc- 
ity depends on y and z coordinates and the meridional 
flow varies only radially) we get the 'dispersion relation' 
(|17|l through introducing an 'effective frequency' because 
of the rather small shear rates and, consequently, a very 
slow and small variation of the wavenumbers and frequen- 
cies in time. The non-modal changes of the 'effective fre- 
quency' and wavevector result in a slight deformation of 
the ray path. In particular, the radial component of the 
wavevector k z now varies along the ray path because of 
two reasons. The first one is the usual variation due to 
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the temperature gradient. The second reason is related 
to the non-modal variation (see Eq. $121 ). If k z increases 
(decreases) within a given layer, when the mode propa- 
gates from the upper turning point to the lower one, then 
the non-modal variation of the radial wavenumber turns 
and extends (shortens) the wavevector. The wavevector 
behaves oppositely when the mode propagates from the 
lower turning point to the upper one: within the same 
layer it turns in the opposite direction and its length re- 
spectively decreases (increases) resulting in a compensa- 
tion of the residual of the 'effective frequency' produced 
within the same layer before. Therefore, the combined ac- 
tion of the refraction and non-modal effects produces the 
shape of the ray path in the kho — k Z Q plane schematically 
shown in Figure ITU1 (panel B). 

It should be noticed that the uncompensated finite 
residuals in modal wavelength (frequency) can arise only 
in the case of nonzero a and c coefficients. This is the 
reason why one can claim that the slight non-modal evo- 
lution of the p-modes occurs only under the joint effect of 
the gradients of the rotational (parameter a) and merid- 
ional (parameter c) flows. The nonzero value a leads to 
the variation of the k y (Eq. (JTTJ) in time. Because of this 
the ray paths do not stay in the kho — k z o plane, but in- 
stead they become 3D curves. The schematic view of these 
curves are shown in Figure ITUI (panel C). 

The residuals occurring during two neighboring reflec- 
tions are very small and they can not be detected in the 
observational data of the modal travel times and distances. 
The non-modal variation of the modal properties is a cu- 
mulative effect in the case of the trapped modes occurring 
only after a significant number of reflections, while the 
effect related to the Doppler shift of the frequencies is im- 
mediately detectable as it systematically shifts all frequen- 
cies at once. Formally speaking the non-modal variation of 
the frequency is a 'second order' effect and it only becomes 
significant after some time after the modal excitation. The 
effect considered in this paper can not be detected by ring 
diagram methods either. As it is well known, the acoustic 
rings represent slices of the 3D local power spectrum at a 
given fixed frequency and deformation of the ring demon- 
strates the shift of frequencies due to advection (for ex- 
ample see Hill (198S|). While because of the non-modal 
variation the 'effective frequency' slowly migrates along 
and/or across the corresponding ridge. That is why we do 
not expect that the non-modal variation will cause any 
additional deformation of the slices of the spectrum by 
fixed frequency planes. However, some signature of the 
non-modal variation can arise because of the azimuthal 
variation of the height (or width of the ridges). But this 
issue could become the subject of a separate study in order 
to verify whether the order of magnitude of the mentioned 
azimuthal variation is compatible with the resolution of 
the employed method (in this case a FFT). To detect the 
non-modal effects efficiently it is convenient to develop a 
different kind of multidimensional spectral method involv- 
ing different methods of data analysis with variable tem- 
poral and spatial scales (such as a window fourier trans- 



form or a wavelet analysis). This latter issue is out of the 
scope of current work and could be a matter of future 
studies. 

6. Conclusions 

The influence of the observed subsurface flow inhomogene- 
ity on solar p-modes has been studied in the framework 
of a non-modal analysis. Particular attention was given 
to the possible role of the non-modal time-dependent ef- 
fects in the formation of the observed p-mode power spec- 
tra. We suggest the hypothesis that several properties 
of the high-degree p-modes could be attributed to the 
slight ('non-modal') deviation of the temporal behavior 
of this kind of modes from a purely ('modal') exponen- 
tial (~ exp(iuit)) evolution. This effect is related with 
the non self-adjointness of the governing equations which 
leads to the existence of alternative solution of these equa- 
tions describing the temporal evolution (in general non- 
exponential) of perturbations. This statement is based on 
the fact that high degree p-modes are mostly confined in 
the upper thin subsurface shell with strongly pronounced 
inhomogeneous flows. We examined this hypothesis and, 
in particular, it was found that: 

a) It is a well-established observational fact that the line 
widths of the peaks in the p-modes power spectrum 
increase with angular degree of modes. In the previous 
section we have discussed the observational aspects of 
the flow inhomogeneity related phenomena. Our re- 
sults show that the nonuniform flows can contribute 
to the excess line widths (compared to the inverse 
observational time) in the observed power spectrum. 
Clearly, the effect addressed in this paper acts along 
with other damping mechanisms. As the frequencies 
in the spectrum are concentrated in the peaks around 
some central frequencies the temporal rates of stochas- 
tic excitement of the modes by different convective 
sources are greater than the characteristic time of the 
significant changes in the modal properties. That is 
why the effect discussed here can not be observed as 
a shift of the peak in the power spectrum. Although, 
some shift can appear because of the different values 
of residuals in frequency of modes propagating in dif- 
ferent horizontal directions with respect to the merid- 
ional flow, but that effect would be much less signifi- 
cant than the changes of the particular wave frequency 
excited by a given convective source. Another conse- 
quence supporting this scenario is that the residuals 
appear due to the asymmetry of the considered veloc- 
ity field. In the case when the terms in the approximate 
'dispersion relation' containing information about the 
meridional flow are negligibly small, the modes un- 
dergo a much smaller net effect from the flow non- 
uniformity during the propagation Sun in the inward 
and outward directions. 

b) We considered the effect of the inhomogeneous back- 
ground flows on the lower degree {t < 200) modes. 
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These modes have their turning points deeper in the 
convection zone, where we have only latitudinal dif- 
ferential rotation. In this case, the effect of inhomoge- 
neous flows is of the order of 6-12 /iHz and this fact is 
in agreement with (and it confirms) the observational 
evidence that the theory of mode excitement / damping 
explains the observed energetic spectrum very well for 
values of the angular degree £ < 200. 
c) Finally, the joint effect of the differential rotation and 
the inhomogeneous meridional flow on the oscillation 
modes with very high angular degree (£ » 200) causes 
significant changes of the mode frequencies in time. 
This may mean that the evolution of the p-modes 
due to the background velocity shear takes part in 
the distribution of the wave energy over frequencies 
and angular degrees. Together with several other pos- 
sible damping mechanisms, this effect can explain or 
contribute to the observed energetic spectrum of the 
modes that shows a systematic deficit of the contri- 
butions from high degree modes. This may be caused 
by partial mixing of power related to the given mode 
with the acoustic noise. Hence, the suggested mech- 
anism may be the cause of the discrepancies between 
theoretical predictions from the model of turbulent ex- 
itation/damping of oscillations on the sun on the one 
hand, and hclioscismic observations of very high degree 
p-modes on the other hand. Additionally, the changes 
of modal frequency due to the non-modal effects might 
be a possible reason for the increases of the wide line 
widths of very high degree solar p-modes. 

The assumptions we made during the analysis of our 
theoretical results enable us to study the main features of 
the modes propagating in the nonuniform flow. But they 
only allow us to understand the 'principle' of the process 
and to estimate the possible changes of the modal proper- 
ties related to the non-modal effects. In order to study this 
problem quantitatively more systematically the following 
efforts have to be made in the future: (i) The analytical 
model should be further developed and probably multi- 
dimensional spectra should be constructed (analogous to 
the ring diagrams) involving also different methods of data 
analysis (say wavelets); (ii) Direct numerical simulations 
of the governing equations in real space and time should be 
carried out; (iii) The results of these simulations should be 
analyzed and compared with the new observational data. 
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Fig. 8. The variation of the modal frequencies in time calculated in the case of the velocity shear profile described 
in Section 4.1 (see Table 1). The direction of the horizontal wavevector is taken as (3 = ±45°. Curves a\ and ai 
(dashed lines) correspond to the sample mode with £ = 95, n = 4 and with basic modal frequency vq = 2.3761 mHz. 
Similarly the 61-62 (dashed dotted lines) and c\-ci (dotted lines) respectively are curves for the modes t = 263, n = 6, 
vq = 4.2336 mHz and £ = 673, n — 2, v$ — 3.8202 mHz. And the curves di-cfo (solid lines) correspond to the same 
mode with I = 951 considered in section 3.2. Different panels show results of calculations in the three different cases 
of the latitudinal gradient of the rotational velocity: Panel A: the parameter a = — 1 • 10~ 8 s _1 ; panel B: a = —1 • 10~ 7 
s" 1 ; Panel C: a = -3 ■ 1CT 7 s" 1 ; Panel D: a = -6 • 1(T 7 s" 1 ; Panel E: a = -1 • 1(T 6 s" 1 . Note that the ordinate 
axes in panels A-E have different scales. In panel F the frequency residuals arising after a 10 hour period are plotted 
as functions of the shear parameter |o| absolute values. 
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Fig. 9. The schematic view of the peak(s) in the power spectrum corresponding to the oscillation mode with a given 
angular degree, a) In the static medium b) in the case of uniform flow; c) Case of the inhomogeneous flow without the 
non-modal effects; d) Case of the inhomogeneous flow including the non-modal effects. 




Fig. 10. The schematic view of the modal ray paths: A) without non-modal effects; B)Same as on panel A (solid 
line) and the non-modal deformation of the ray path due to the temporal evolution of the radial component k z of 
the wavevector, which first increases and then decreases in length (dashed line) and vice versa (dash-dotted lines); 
C) The non-modal deformation of the ray path because of the temporal evolution of the meridional component k y of 
the wavevector (dashed and dash-dotted lines). The modes do not stay in one plane, but the trajectories become 3D 
curves. 



